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ABSTRACT 

Context. We present a multi- wavelength, high-resolution observational survey of the young binary system Haro 6-10 (GV Tau, IRAS 
04263+2426), which is harbouring one of the few known infrared companions. 

Aims. The primary goal of this project is to determine the physical and geometrical properties of the circumstellar and circumbinary 
material in the Haro 6-10 system. 

Methods. High-resolution optical (HST/WFPC2) and near-infrared (VLT/NACO) images in different bands were analysed to investi- 
gate the large-scale structures of the material around the binary. Mid-infrared interferometry (VLTl/MlDl) and spectroscopy (T1MM12 
at the 3.6m ESO telescope) were carried out to determine the structure and optical depth of the circumstellar material around the in- 
dividual components. 

Results. The multi- wavelength observations suggest that both components of the binary system Haro 6-10 are embedded in a com- 
mon envelope. The measured extinction indicates a dust composition of the envelope similar to that of the interstellar medium. Each 
component of the system has a circumstellar disc-like structure typical of young stars. The discs are highly misaligned: the northern 
component is seen almost edge-on and the southern component is an almost face-on disc. 

Conclusions. The two main formation scenarios of binary systems with misaligned discs are the gravitational capture of a passing 
object in a dense environment, and the fragmentation of the collapsing molecular cloud. Given the low-density environment of the 
Taurus-Aurigae star-forming region, the first scenario is unlikely for Haro 6-10. The binary system most probably formed via frag- 
mentation of two different parts of the collapsing molecular cloud combined with other dynamical processes related to the cloud 
and/or the protostars. This can be the explanation also for other binary systems with an infrared companion. 
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(«— ») 1. Introduction ularcloud, this is not the case for rotational fragmentation, which 

tH- is caused by asymmetric instabilities in a rotating disc or ring. If 

^^ The formation of binary and multiple systems seems_to_be jj^^ proto-binary accretes (only) gas with low angular momen- 

a^mmon phenomenon in star formation (e.g., MathieuelalJ ^^^^ ^ ^-^^ -^ ^^^^^ ^^ly ^^^^^^ ^1^^ p^.^^^^^^ ^^^^ ^l^il^ ^^ 

r- ; 12000). Surveys of young star-forming regions with high spatial ^g^ondary may have none. If instead the seco ndary developed 

2 . resolution allow a statistical description of the binary parame- -^^ ^^^ ^-^^^ ^j^^ tj^^ j^ ^^ j^^^^ ^ ^-^^ (iBate & Bonnelll 

'^ ters and a comparison betwe en regions of yariotis ages and with IT9Q7K 
". . different environments (e.g.. iLeinert et alj|1993l : iDuchene et alj 

> 112007). Binarity has been studied also in young stellar clusters According to the two different scenarios of binary forma- 

k> ' of various ages to inve stigate its relation to the evolution of pro- tio"' the resulting circumstellar discs may have different incli- 

/\ ' toplan etary discs (e.g.. lCIarke & PringIell991btlBouwman et"d] nations. When the binary is formed via rotational fragmenta- 

^4 ' 120061) tion the discs are expected to be preferentially aligned. If, how- 

■ ■ ■ ' Different scenarios have been proposed to explain the formation ever, the stellar components and their own discs form direcdy 

of binai-y systems. The standai'd scenario supports the hypoth- from a collapsing core ('prompt initial fragmentation' : Pringld 

esis that a binary or multiple system forms when the core of a 1989), this can lead to two misaligned discs (jBate et al.||2000|). 

molecular cloud fragments during its gravitational collapse. The Under more complicated initial conditions even a small clus- 

fragmentati on process can be divided into two main classes : di- ter can be formed from such a core. Moreover, interactions in 

rect (e.g.. lBoss & Bodenheim er 1979; Bate & Burkert 1997) and these clusters between stars with discs can cause a dissipative 

rotatio na l fragmentation (e.g ., Bonnell .1994: Bonnell & Bate encounter (e.g., Clarke & Pringle 1991b). This may lead to the 

[T994allbl; iBurkert et all [1991 . While the direct fragmentation formation of binary or mukiple systems via capture and thus to 

strongly depends on the initial density distribution in the molec- binaries with misaUgned discs. In dense star-forming regions, a 

mas sive accretion disc may e ven cause the capture of a passing 

Send offprint requests to: V. Roccatagliata: rocca@stsci.edu stai' (jCIarke & Pringleil991a|). On the other hand, strong tidal in- 

* Based on observations collected at the European Southern teractions can lead to aligned discs. Another possibility to have 

Observatory, Paranal, Chile (Proposal IDs: 070.C-0544(A); 072.C- misaligned discs is that the gravitational interaction of a passing 

032 1(A); 274.C-5032(A); 076.C-0200(A)) object causes the disc's tilt (and hence also the direction of a pos- 
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sible jet) in a binary system dBate et alJl2000l) . lHelle3(ll993h es- 
timated that in one million years, in a small cluster, a few percent 
of the stellar systems can experience a tilt of a few degrees dur- 
ing the star encounters. The effect of stellar encount ers in differ- 
ent clu ster environments has been investigated by Olczak et alj 
(I2OIOI) . They found that the encounter rate remains unaffected by 
the size of the stellar population, while it strongly depends on the 
cluster density: in clusters less dense and less massive than the 
ONC, massive stars dominate the encounter-induced disc-mass 
loss, whereas in denser and more massive clusters the disc-mass 
removal is controlled by low-mass stars. 

In systems near the end of their accretion phase, the infall of 
a small amount of material with a different angular momentum 
to th e orbit can cause a misalignment of discs originally aligned 
(e.g.. lBateetani2000l) . 

Therefore, from the alignment of circumstellar discs it is pos- 
sible to understand the dynamics that most likely took place dur- 
ing the star-formation process. 

Studying disc alignment requires observations of proto- 
planetary discs around young binary stars with high spa- 
tial resolution. Previous st udies deterniined t he inclination b^ 
means of po larimetry (e.g. jMonin et al.lfT998H Wolf et al. 2001 



Table 1. Log of TIMMI 2 observations of Haro 6-10 and of the 
spectro-photometric standard stars. 



Monin 



T99a 



etal.ll2006>) 



orie ntation of the jets (e.g., .Eisloffel et al 
Davis etal.1 II 997h. and direct observati ons of the discs 
(e.g., 'Koreskol Il998t IStapelfeldt et al.l [1998b'). Another tech- 
nique, that allows to directly measure size, inclination and orien- 
tation of the inner parts of protoplanetary discs is long-baseline 
interferometry at infrared wavelengths. A recent interferometric 
study of T Tau sho wed that the three discs in this triple sys- 
tem are misaligned (Ratzkaet al. 2009). Therefore, the forma- 
tion process of T Tau probably has been highly dynamic. 

We present in this paper a multi-wavelength, high-resolution 
observational survey of the young binary system Haro 6-10 
(GV Tau, IRAS 04263+2426). Haro 6-10 is a T Tauri binary sys- 
tem residing within the L-1524 molecular cloud. It is composed 
of an optically visible source, Haro 6-10 S, and an infrared com- 
panion (IRC), Haro 6-10 N. With an angular separation of l."2, 
Haro 6-10 N is actually located -160 AU north of Haro 6-10 S at 
the distance of the Taurus- Auruga star- forming region. The bi- 
nary system was firstly resolved bv lLeinert & HaasI (11989 ) using 
speckle interferometry. Haro 6-10 N has a very red spectral en- 
ergy distribution and is brighter than Haro 6-10 S at wavelengths 
longer than 4 fim. The binary system is variable, particularly in 
the n e ar-lR, on a timesca le as short as a month (Leinert et aL 
l200lh . lLeinertetal.l (l200lh suggest that the variabiHty of Haro 6- 
10 S is associated to inhomogeneities in the circumstellar mate- 
rial, while the variable behaviour of the northern component is 
likely associated to a variable accretion mechanism. The first 
resolved images in t he near-infrared of this system, obtained by 
iMenard et al.l(ll993h . suggest the presence of a flat circumbinary 
envelope, or disc. 

The infrared spect rum of Haro 6-10 b etween 8 and 13 fim 
has been interpreted bv lHanner et al.l (Il998l) by an optically thick 
plus an optically thin emission, where the optically thin emis- 
sion could be associated to an envelope, while the optically thick 
emission is associated to a disc, to the star itself, or to a fea- 
tureless dust continuum. The Spitzer/IRS spectrum (between 5- 
36 fxm; Furlan et al. 2008) shows a broad silicate absorption fea- 
ture at 10 fim, but the spatial resolution of the instrument is too 
lo w to resolve the binary. 

iReipurth et al.l (l2004l) obtained high-resolution 3.6 cm ob- 
servations with the VLA. In the 3.6 cm map the two compo- 
nents of the binary are spatially resolved. At this wavelength 
the radio continuum emission originates from shock-ionised gas 



Date 



UT 



Objects Airmass 



2002 Feb. 06 


00:10 


HD32887 


1.02 




00:41 


Haro 6-10 


1.71 




00:53 


Haro 6-10 


1.73 




01:41 


HD29139 


1.56 




07:46 


HD123139 


1.11 


2002 Dec. 26 


00:42 


HD29139 


1.82 




01:35 


Haro 6-10 


1.84 




01:51 


Haro 6-10 


1.78 




05:00 


HD32887 


1.07 



very cl ose to the out fl ow ex citing source (e.g., 'Angla da et al.l 
119981) . IReipurth et aP (12004 ) suggest that Haro 6-10 S is the 
driving source of the large outflow activity seen in the near- 
infrared. HCN and C2 H2 in absorption to wards Haro 6-10 N 
were firstly detected by iGibb et alJ (l2007l) in the near-infrared 
using NIRSPEC at the Keck Telescope and one year later by 
iDoppmann et al.l (l2008h . using the same instrument. 

In this paper, we are using a multi-wavelength approach, in- 
cluding different techniques, to characterise the binary system 
Haro 6-10. This paper is organised as follows: in Section|2]we 
present the archival optical and near-infrared photometric obser- 
vations, and our mid-infrared spectroscopic and interferometric 
observations. The data reduction and the results are described in 
Section 13 The analysis and the discussion of the results are pre- 
sented in Sections |4] and |5] respectively. Conclusions are drawn 
in Section|6] 



2. Observations 

In this section we present the observations collected from the 
archives and our new TIMMI2 and VLT/MIDI observations of 
Haro 6-10 N and S. 



2.1. Archival Optical and Near-infrared images 

HST imaging observations of Haro 6-10 were obtained as part 
of the WFPC2 snapshot survey of nearby T Tauri stars (HST 
7387 program, PI: K. Stapelfeldt) in November 1998, using 
the HST PCI (Planetary Camera 1) and the HST WF3 (Wide- 
Field Camera 3), which have a spatial scale of 0.05"px ' and 
O.r'px"', respectively. The observations have been carried out 
with the F814W and F606W filters and consist of one expo- 
sure of 120 seconds with the F814W filter, and 2 exposures with 
the F606W filter of 200 seconds and 600 seconds, respectively. 
The calibrated data are available from the ESO/ST-ECF Science 
Archive FacilitjQ. 

The archival near-infrared NACO/VLT observations of 
Haro 6-10 N and S, obtained in December 2003 (Program 
ID: 072.C-0321(A), PI: J.Bouvier), have been also included in 
our analysis. These data are also available from the ESO/ST- 
ECF Science Archive Facility'. While the HST images are un- 
published, the binary parameters derived from the NACO ob- 
servations were already published, but only in the context of 
an high-resolution survey of multiple young stellar systems 
(Duchene et al. 2007) . 



http://archive.eso.org 
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Table 2. Log of MIDI observations of Haro 6-10 N and S and calibrators observed, 
visibilities and fluxes; **: quality not good enough to be analysed 



calibrators used to obtain the final calibrated 



Airmass 



Date 



UT 



Objects 



Proj. Baseline 
[m] 



Pos. Angle 
[deg] 



Notes 









[UT2-UT3] 






2004 Dec. 28 00:18 


HD37160 


2.08 


27.61 


14.61 


Calibrator* 


01:01 


Haro 6-10 N 


1.73 


28.43 


48.67 


Science 


02:38 


Haro 6-10 S 


1.60 


38.13 


52.15 


Science 


03:10 


HD37160 


1.22 


39.98 


43.62 


Calibrator* 


04:08 


HD37160 


1.21 


43.61 


46.24 


Calibrator* 


05:26 


HD50778 


1.03 


46.43 


42.72 


Calibrator 


06:10 


HD94510 


1.43 


45.17 


15.53 


Calibrator 


07:02 


HD49161 


1.40 


46.61 


45.68 


Calibrator 


07:45 


HD94510 


1.25 


44.00 
[UT3-UT4] 


32.44 


Calibrator 


2005 Nov. 23 04:17 


Haro 6-10 N** 


1.56 


62.13 


106 


Science 


05:22 


Haro 6-10 S 


1.54 


58.60 


95.6 


Science 


05:41 


Haro 6-10 N 


1.56 


56.80 


92.6 


Science 


06:09 


HD25604 


1.61 


50.53 


85.4 


Calibrator* 


07:49 


HD33042 


1.23 


61.60 


116.5 


Calibrator 
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Fig. 1. Perpendicular (to the dispersion direction) cut of a 
TIMMI 2 bi-dimensional spectrum obtained on December 26 
along with the double Gaussian function fit utilised for the spec- 
tra extraction. Data are shown as solid line, fit with diamonds. 
The fit to each component is shown as dotted line for the north- 
ern component and dashed line for the southern component. The 
TIMMI2 pixel-scale in the N-band was 0.45". 



2.2. Mid-infrared spectroscopic and interferometric 
observations 

Haro 6-10 was spectroscopically observed in the mid-infrared 
with the TIMMI2 instrument mounted at the 3.6 m ESO tele- 
scope at La Silla (Chile) during two campaigns in February and 
December 2002, see Table [1] The TIMMI2 pixel-scale in the N- 
band was 0.45". The slit width used during the observations was 
3.0"and 1.2"in February and December, respectively. 

We obtained two mid-infrared interferometric measurements 
of both Haro 6-10 S and Haro 6-10 N with MIDI at the VLTI 
during the Guaranteed Time Observations (GTO) in December 
2004 and November 2005. The first observations were carried 
out using the baseline UT2-UT3, and the second observations 
with the baseline UT3-UT4. Details of the observations are listed 



in Table 12] The MIDI slit width was 0.5". The interferometric 
observations with MIDI were performed using the standard pro- 
cedure for the high sensitivity mode which consists of 1 .) acqui- 
sition images taken in the lower N-band around 8.7 fj.m, 2.) mea- 
surements of the spectrally dispersed correlated flux, 3.) single- 
dish spectra from each telescope. Spectrally resolved visibilities 
are then derived as the ratio of the correlated flux and the single- 
dish spectra. The spectral resolution obtained in the observations 
of Haro 6-10 is A/AA ~30 as the prism was used as dispersive 
element. A detaile d description of the M IDI observing sequence 
is reported, e.g., in lRatzka et al.l (l2009i) . 

Soon after each scientific observation, the sequence is re- 
peated for a standard star of known diameter These calibra- 
tors are required to determine instrumental and atmospheric ef- 
fects. The calibrators, partly spectro-photometric, also served to 
calibrate the absolute flux as well as the correlated fluxes. We 
checked the stability of the transfer function using all the calibra- 
tors observed during each night. The final calibration has been 
done using the calibrators closer in time to the science target (see 
TablelUl. 



3. Data reduction and results 

3. 1 . Optical HST pliotometry 

We used the calibrated HST data of Haro 6-10 already processed 
by the standard data reduction pipelinfl which are available at 
the ESO/ST-ECF Science Archive Facility. Cosmic rays in each 
single exposure hav e been removed using a specific IDL code 
(Ivan Dokkumll200ll) . Haro 6-10 S is saturated in the long expo- 
sure of 600 sec in the F606W filter. Figure[3]shows a stack of the 
F606W 200s and 600s exposures (with Haro 6-10 S replaced by 
the non-saturated pixel in the 200 sec exposure), and the 120 sec 
exposure obtained with the F814W filter 

Aperture photometry of the short exposure in both filters 
(F606W and F814W) was carried out on Haro 6-10 N and 
Haro 6-10 S. An aperture of 0.50"(corresponding to 5 pixels in 
the F606W filter, 10 pixels in the F814W filter, respectively) has 
been adopted for the aperture photometry of the bright object 



http://archive.eso.org/archive/hst/ 
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Table 3. Calibrated magnitudes obtained with HST and 
VLT/NACO. 





Haro 6-10 N 


Haro 6-10 S 


flux 




[mag] 


[mag] 


ratio 


V (F606W) 


23.445±0.084 


17.169±0.003 


0.003 


I (F814W) 


2 1.041 ±0.059 


15.012±0.002 


0.004 


H 


14.2+0.1 


11.1 ±0.1 


0.06 


Ks 


10.7 ±0.1 


9.2 ±0.1 


0.24 



Haro 6-10 S. The sky was estimated in the two filters in an annu- 
lus between 0.55"and 1.55"and 0.6"and 1.2", respectively. An 
aperture of 0.2"has been adopted for the aperture photometry of 
the faintinfrared companion Haro 6-lONin both filters. Due to a 
non-perfect cleaning of the cosmic rays and due to saturation of 
the central pixel, the long 600 sec exposure in the F606W filter 
cannot be used for aperture photometry. 

The final magnitudes have been computed, including the 
gain ratio coiTection of the PCI an d WF3 cameras and using 
the zero points and gain ratios from lHoltzman et aP (Il995h . The 
values obtained are listed in Table [3] 



3.2. Near-infrared photometry with NACO. 

NACO observations of Haro 6-10 were reduced with the data 
reduction pipeline provided by ESCQ. Both components of the 
binary system were saturated in the L-band exposure. We de- 
rived the zero-points of the photometry from standard stars ob- 
served at the beginning and the end of the night. An aperture of 
0.2" (corresponding to 15 pixels) has been adopted for the aper- 
ture photometry and the sky was estimated in an annulus be- 
tween 0.30"and 0.35". The results of the aperture photometry in 
the H- and Ks-bands are reported in Table [3] 



3.3. TIMMI2 mid-infrared spectroscopy 

The northern, infrared component brightens at 1 3 yum from about 
20 Jy in February to about 30 Jy in December 2002. The southern 
component tends to become fainter at short wavelengths with a 
less pronounced or maybe broader silicate absorption feature. 
The derived spectra show significant temporal changes in the 
overall shape and absolute flux levels. Individual spe ctra of both 
binar y components have been already presented in iPrzygodda 
(1200 4). The spectra of Haro 6-10 S and its infrared companion, 
Haro 6-10 N, are only separated by a few pixels on the TIMMI 2 
detector: For this reason we reanalysed the data of February 6 
and December 26 with a different method to check how large the 
effect of the extraction method on the spectral variations is. With 
the aim of optimising the spectrum extraction we have first fitted 
the preprocessed bi-dimensional TIMMI 2 spectrum with a dou- 
ble Gaussian for each wavelength column (Figure[T]l. In this way 
the two spectra of each binary component are traced. Besides the 
fitted spectra, the fitting routine delivers synthetic images that 
can be compared to the original files. We found that, except the 
wings of the lines, the results of the fits are of good quality. The 
calibrators were also fitted and extracted with a Gaussian func- 
tion to avoid any systematic error. 

The correction for telluric absorption as well as the abso- 
lute flux calibration was achieved by means of observations 
of spectro-photometric standard stars observed at different air- 



masses during each night; the extracted spectrum of each stan- 
dard star was flux-calibrated using the spectral templates by 
[Cohen et al. ( 1999). Finafly the spectra of Haro 6-10 S/N were 
scaled with the flux calibrated standa rd star spectra. Th e wave- 
length dispersion solution derived bv Przygoddal (|2004|) was ap- 
plied. The final TIMMI 2 spectra are shown in Fig.|2] 

The spectra of the southern source are found to be stable 
in the course of 2002. This ensures the reliability of the new 
extraction method used. The TIMMI2 spectra show a similar 
shape than the MIDI spectrum taken two years later (details in 
section 13.41 1. but the absolute level is about three times larger. 
The TIMMI2 spectra of the northern component sh ow changes 
in flux in the upper part of the N-band (also shown by iPrzvgoddal 
120041) . However, the MIDI spectrum of Haro 6-10 N resembles 
the TIMMI2 data from February 2002 very well. The instru- 
mental FWHM increases from 0.46"at 8 yum to 0.74"at 13 yum. 
The effective FWHM, measured directly from the spectrum is 
0.9"and does not depend on the wavelength. We notice that 
the differences between the TIMMI2 and MIDI spectra con- 
cern only the southern component of the binary system, Haro 6- 
10 S, which is the fainter one. Although the different FWHM 
of the two instruments, TIMMI2 and MIDI, might introduce in 
the TIMMI2 spectra background emission from the suiTound- 
ing nebulosity, this would have affected both the spectra of Haro 
6-10 N and S; however, since the northern component is about 
10 times brighter than the southern one, the contamination from 
the background might significantly increase only the flux of the 
fainter component Haro 6-10 S. Another problem might be that 
the deblending of the two TIMMI2 spectra is is highly contam- 
inated by the northern brighter source. However, the latter can 
probably not alone explain the large discrepancy between the 
MIDI and the TIMMI2 spectrum of Haro 6-10 S. 

3.4. MIDI Data 

The data reduction of interferometric data with MIDI has been 
performed with the MIAh-EWS 1.5 data reduction software pro- 
vided by the MIDI consortiurrQ The results of the data reduction 
are the raw correlated spectrum Fj:™ and the total (single-dish) 
raw spectrum Fiaw We here use the results of the EWS package, 
which determines the correlated flux FJ:™ based on a "coherent" 
algorithm. MIA, based on power spectrum analysis, was used to 
confirm the results. 



3.4.1. Visibility 

We compute the raw visibihties from 



'^raw — 



rcoiT 
raw 



(1) 



where Fraw is the geometrical mean of the total fluxes of 
the incoming beams from telescope A and B, i.e., F^v, = 

V ^ raw ^ raw ■ 

The calibrated visibilities Vca\ of the science object are ob- 
tained by dividing the raw visibilities, Vi-aw, by the instrumental 
visibilities derived from one or more calibrators observed within 
the same night, usually the closest in time. In Table|2]we compile 
the calibrators observed during each night, marking those, ob- 
served closest in time to the science target, that have been used to 
calibrate the final visibilities. The other calibrators were used to 



http://www.eso.org/sci/data-processing/software/pipelines/naco/ 



"^ http://www.strw.leidenuniv.nl/~koehler/MIA+EWS- 
Manual/index.html 
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Fig. 2. TIMMI 2 spectrum of Haro 6-10 N (IRC) (top) and 
Haro 6-10 S (below). The grey highlighted pattern correspond 
to strong atmospheric O3 and CO2 absorption bands at 9.6 and 
12.5 /zm, respectively. The measurements in these regions are 
not considered in our following analysis. 



check the stabiUty of the transfer function. The error associated 
to the calibrated visibility of the science targets is the standard 
deviation resulting from the calibration with different calibrator 
stars for each night. 



3.4.2. Spectroscopy 

As part of the observing scheme, MIDI acquires the N-band 
spectrum of the target (spectral resolution equal to that of the in- 
terferometric signal) with both telescopes separately. MIDI spec- 
tra were extracted using the MIA-i-EWS software. In this regards 
we used the three observations of the spectro-photometric stan- 
dard star HD 37160 obtained in December 2004 (see Table |2]i. 
In Fig.|2]we compare the MIDI and TIMMI2 spectroscopy. The 
results are very similar for Haro 6-10 N, while for Haro 6-10 S 
TIMMI2 yields a substantially higher absolute flux. This is likely 
due to the strong effect of the deblending on the fainter compo- 
nent and/or to the higher spatial resolution of the VLT compared 
to the 3.6m telescope. The TIMMI2 spectrum of Haro 6-10 S 
might be contaminated by the northern component (see also dis- 
cussion in Section [33] ) or by large-scale background emission. 



3.4.3. Correlated flux 

To flux calibrate the correlated spectrum one approach is to use 
the following equation 



cal 



Vc 



cal 



call 



(2) 



where Vcai is the calibrated visibility and Fcai is the flux- 
calibrated spectrum. An alternative approach is to directly cal- 
ibrate the raw correlated flux, Fj:™, by means of spectro- 
photometric standards. The results obtained using the two meth- 
ods for bright sources agree well, while there are some discrep- 
ancies for the faint source Haro 6-10 S. The correlated fluxes for 
our analysis are all consistently computed with Equation (|2]l. 

4. Analysis 

In this section, we present the analysis of our multi-wavelength 
observations of Haro 6-10. We analyse first the large-scale struc- 
ture of the binary by investigating visual and near-infrared im- 
ages as well as mid-infrared spectra. In a second step, we exam- 
ine the geometrical properties of the protoplanetary discs around 
the binary components through modelling MIDI mid-infrared 
visibilities and correlated fluxes. 



4.1. Large-scale morphology in the optical and near infrared 

Figure |3] shows the WFPC2/HST optical images and the 
VLT/NACO near-infrared images. Both components of the bi- 
nary system are resolved from the optical to the mid-infrared 
wavelengths. While in the optical the infrared companion is only 
barely detected, it becomes brighter in the near-infrared. The flux 
of Haro 6-10 N is 0.3% of the flux of Haro 6-10 S in die visual 
and increases up to 24% in the K band. The photometric results 
and the flux ratios are listed in Table [3] In the N-band Haro 6- 
10 N is brighter than the southern companion (Figure|2|i. 

The resolution of 0.05"/px and 0.1 "/px of the HST images 
allows to disentangle structures and the filaments in the en- 
vironment of the Haro 6-10 binary system. An arc-like struc- 
ture is clearly detected up to a distance of ~6"S and ~3"W 
from the southern component, where it fades, together with 
other north-south oriented filamentary struct ures, in the underly- 
ing nebulosity. Pre vious observations (e.g.. iLeinert et al.l 120011 : 
iKo resko et al.l 11999 ) discovered the short and bright arc-like 
structure east of Haro 6-10 S (see Section ISTt . Another promi- 
nent structure is a bright knot less than l"south-west of Haro 6- 
10 S. In the H- and K-band these filament structures become 
fainter and are undetectable at longer wavelengths. The southern 
component appears to be barely resolved in the E-W direction. 

4.2. Optical extinction and dust composition 

To derive the optical depth of the circumstellar material we fit- 
ted the continuum with a first-order polynomial to the TIMMI2 
spectrum between 8 and 13 jum, assuming that at these wave- 
lengths only the continuum contributes to the observed spec- 
trum. The depth of the absorption feature was then determined 
with a secon d-order of polynom ial fit. A similar approach was 
presented bv lOuanz et al.l (l2007b . We compute the optical depth 
^"silicate using the following equation: 



And from this we derive the extinction AsHicate by 
Asiiicate = -2.5 ■ logioie'^"""'°) ■ 



(3) 



(4) 
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Fig. 4. Normalised spectra of Haro 6- 1 N (top) and Haro 6- 1 S (bottom) obtained with TIMMI2 on February and December 2002 
and with MIDI on December 2004, with the computed maximum optical depth. A second-order polynomial has been fitted to the 
data between 8.8 and 10.2 jum (dot-dashed lines) to compute the maximum optical depth (r). For comparison the ISM silicate profile 
from lKemper et all (|2004|) is overplotted (dashed red lines). 



The optical extinction Ay was computed using the relation 
Av/Tgj - 19.3 from Mathis ( 1998). All the results ai-e listed in 
Table|4] We derive an average visual extinction of 18.7±4.1 mag 
for Haro 6-10 N and 13.0+3.8 mag for Haro 6-10 S. Our val- 
ues arejiigherbut_wifli a similar ratio than the values obtained 
by iLeinert et al.l ( 1200 li) . From the ice-absorption feature they 
derived 10.6+0.8 mag for Haro 6-10 N and 7.0±0.7 mag for 
Haro 6-10 S. For such a young system, this absorption is most 
probably caused by a circumbinary envelope. The spectra of 
Haro 6-10 N and Haro 6-10 S are compared in Fig. |4] with the 
profile of the silicate dust present in the interstellar medium 
(ISM). This profile was measured by Kemper et al. (2004) and 
is shown as dashed line. The spectra of the binary components 
are very similar to the ISM profile. However, the ISM absorption 
looks broader compared to the absorption observed towards the 
binary. The quality of the spectra does not allow to conclude 
whether these deviations are due to a variation in the chemi- 
cal composition or on the grain size of the dust around the bi- 
nary system. A further analysis of the dust composition using, 
e.g., a simple fitting of th e N band spectrum (as done for Herbig 
and TTauri stars, e.g., bv lBouwman et alll200ll: ISchegerer et al.l 
12006) is beyond the goal of this paper. 



4.3. Small-scale morphology: Geometry of the discs 

In this section we introduce the modelling strategy adopted to in- 
terpret the MIDI interferometric measurements. The characteris- 
tic size of the emitting region was firstly computed assuming a 
Gaussian brightness distribution. In a second step, the emitting 
sizes are computed fitting the two observations available per ob- 
ject with a Two Black-Body model with a Gaussian distribution 
of the emitters. An alternative approach would be to do the si- 
multaneous modelling of the spectral energy distribution (SED) 
and the mid-infrared interfe rometric measurements, using a ra- 
diative transfer model (e.g.. ISchegerer et al.ll2008t iRatzka et al] 
I2OO9I: ISchegerer et al.ll2009l) . This has the advantage to repro- 
duce in a more realistic way the physical structure of the circum- 
stellar disc. However, the degeneracy in the parameter space left 
when using only a few visibility measurements can be reduced 
only by taking supplementing information into account, e .g., 
high-resolution images at different wavelengths (e.g.. lWoIf et al] 
2003. 2008; Sauter et al. 2009i) . 



4.3.1. Gaussian brightness distribution 

To estimate the size of the emitting region from the observed 
visibility on the basis of a Gaussian brightness distribution is 
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Table 4. Optical depths and extinction values as derived from 
the spectra in Fig. ID 



Object 


Date 


"^silicate 


■^silicate 


-'T^ silicate 


Av 


Haro 6-10 N 


2002 Feb. 


1.0 


9.7 


1.1 


18.8 




2002 Dec. 


1.2 


9.7 


1.3 


22.8 




2004 Dec. 


0.8 


9.7 


0.8 


14.6 


Haro 6-10 S 


2002 Feb. 


0.5 


9.5 


0.5 


8.8 




2002 Dec. 


0.7 


9.7 


0.8 


14.3 




2004 Dec. 


0.8 


9.3 


0.9 


16.0 



Table 5. Derived emitting size in [AU] using the FWHM of a 
Gaussian brightness distribution. 



[/im] 


Haro 6-10 N 

UT2-UT3 UT3-UT4 

[28 m] [62 m] 


Haro 6-10 S 

UT2-UT3 UT3-UT4 

[38 m] [58 m] 


7.9-9.8 

10.0-12.0 

12.1-13.0 


3.4±0.6 
4.3±0.1 
4.6±0.8 


2.5±0.3 
3.2±0.1 
3.5±0.2 


1.8±0.4 
2.9±0.2 
3.2±0.6 


1.0+0.2 
1.4±0.3 
2.2±0.3 






Fig. 3. Upper panels: Color-composite image of Haro 6-10 in 
logarithmic scale. The blue, green, and red colours coiTespond 
to the photometric filter F606W (HST) and the H and K, band 
(NACO). North is up. East is right. Lower panels: HST images 
obtained using the F606W (~V band) and F814W(~I band) fil- 
ters and NACO images in the H, K,, and L'-band. 



a reasonable first approximation alread y applied for bare ly re- 
solved (high visibility) objects (e.g., by 'O uanz et al.ll2006l) . The 
visibility of such a distribution is computed as 

V(/) = exp(-3.56/02) , (5) 

where is the Full Width Half Maximum (FWHM) of the 
Gaussian distribution in arcsec, and / the spatial frequency in 
arcsec"'. For each baseline the emitting sizes of the two discs 



(using the FWHM) are reported in Table |5] as function of wave- 
length. These values are given in astronomical units (AU) as- 
suming a distance of 140 pc. We average the results obtained 
within 3 wavelengths bins: 9.8 - 7.9 yum, 12.0 - 10.0//m and 13.0 
- 12.1 jum. The associated errors are the standard deviations. We 
note that the emitting size increases with wavelength. Such a 
behaviour is expected for thermal disc emission, where the tem- 
perature of the disc decreases with the distance from the cen- 
tral star. This implies that we are resolving two extended struc- 
tures around the central star. Despite the large appearance dis- 
crepancy of the two sources in the visible, the derived sizes of 
the two dusty components in the mid-IR are very similar, even 
when taking into account the uncertainties on the position angle 
of each disk. The three wavelengths bins reflect the three regimes 
of opacity: the 'hot' continuum, the increase of size due to the 
silicate opacity, and then the 'cool' more extended continuum 
near 13 micron. 

4.3.2. The 'Two Black Body' Model 

We note that the change in the slope of the visibility with wave- 
length may suggest the presence of two different emitting re- 
gions contributing to the mid-infrared flux (see Figure |5]). We 
introduce here a more complex (geometrical) model which takes 
into account also the spectrum and correlated flux of the emit- 
ting source. We model simultaneously the MIDI spectiTim (in- 
tegrated over the beam width of 0.5") and the correlated fluxes 
corresponding to spatial scales at lOjum of 0.07" and 0.03" for 
the 28 m and 62 m projected baselines, respectively. 

In a circumstellar disc, in first approximation, we expect a 
radial decrease of the temperature as a function of the distance 
from the central star. Taking into account the vertical extension 
of the disc, a decrease of the disc temperature is expected also 
from the surface to the mid-plane. The mid-infrared flux is there- 
fore expected to come from the inner part of the disc close to 
the central star. In this simplified view we approximate the emit- 
ting source with two black-bodies of different temperatures: each 
black-body is characterised by a Ga ussian brightness distribu- 
tion (a similar approach was used bv lJaffe et alJl2004l) . In order 
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Table 6. Best fit parameters of the 'two black body' model. 
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Fig. 5. Calibrated visibiUties of Haro 6-10 N (top) and Haro 6- 
10 S (bottom). 



to reproduce the dust absorption feature between 8 and 13 fim 
the silicate absorption profile is included. In particular, we use 
the absorpt ion profile d erived from ISC|j towards the centre of 
the galaxy (iKemper et al...2004.) . which is consistent with the ab- 
sorption profile of Haro 6-10 N and S (see Section \42\ . The 
details of the model used are presented in Appendix IBJ 

We used position angles of the discs (Haro 6-10 N: 140 ", 
Haro 6-10 S: 120°) as d erived from the tw o different jets origi- 
nating from each source (iBecket al.l2010l T Beck, private com- 
munication). This has been done to ensure that the number of 
free parameters is not biasing the results. 

A first best model was found minimising the reduced ;if^ (x^ 
divided by the number of free parameters of the fit). In order to 
have a more reliabl e result, following the approach presented by 
iFedeleet all ( 120081) . 

1 . we created 1000 random new data-sets of measured visibili- 
ties around the observed visibilities, assuming a normal dis- 
tribution of the visibility errors; 

2. we computed the best-fit model for each new data-set. 

The final best fit parameters and associated uncertainties corre- 
spond to the mean and the standard deviation of the 1000 fits. 
Fig. |6]shows the histograms of the best fit parameters computed 
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Haro 6-10 S 


r, [AU] 


1.5±0.5 


1.0+0.5 


r. [AU] 


10±2 


7±3 


T, [K] 


900±100 


900±300 


T2[K] 


150±50 


100±50 


in 


80±10 


10+5 



Infrared Satellite Observatory 



for each fit, and the final best fit values are listed in Tabled 
For Haro 6-10 N we found that the emission within 1.5 AU from 
the central star has a temperature of 900A' and a second colder 
emission at 15QK originates within 10 AU from the central star. 
For Haro 6-10 S we found an emission within 1.0 AU from the 
central star with a temperature of 900K and a second colder 
emission at lOOK, originating within 7 AU from the central star 
The disc of Haro 6-10 N is seen close to edge-on (/ = 80°), while 
Haro 6-10 S is almost seen face-on (/ = 10°). The derived incli- 
nations of the binary components are consistent with the fact 
that the central star of Haro 6-10 S is visible in the optical, while 
the star of Haro 6-10 N is obscured by an additional extinction 
represented by the circumstellar disc and is only barely visible 
in the optical. These results are also consistent with the general 
appearance of the binary system described in Section |4jTJ 

In the case of Haro 6-10 N, the model fits the observations 
very well, while for Haro 6-10 S the model fitting is less accu- 
rate; the model reproduces the mid-infrared spectrum, the corre- 
lated spectra, and the absolute value of the visibilities, but does 
not fit their shape. We thus additionally tested the robustness 
of our results to be sure that the Monte Carlo simulations did 
not only find a local minimum in the parameter space instead 
of the absolute minimum. Since the best defined parameters are 
the disc inclinations, we changed their initial values, for Haro 6- 
10 N to a face-on disc (/ = 3°) and for Haro 6-10 S to an edge-on 
disc (/ = 80°). Also in this case the Monte Carlo simulations con- 
verge to the results previously obtained. We also run the whole 
modelling with the position angles as additional free parameters. 
Although the derived position angles differ from the observed 
ones by ~ 70" for Haro 6-10 N and ~ 50° for Haro 6-10 S, the 
results for the remaining parameters are almost identical to the 
above given values. 

Given the higher uncertainties of the measurements for the 
southern component and the poorer fits (Fig. [8]l, the formal er- 
rors derived for its disc parameters are most probably a lower 
limit. The simple size estimates derived in Section l4.3.1l are also 
consistent with higher inclinations (Table|5]l. We further want to 
emphasise, that our approach strongly relies on the assumption 
that the components of Haro 6- 1 can be described by a superpo- 
sition of two Gaussian brightness distributions. If the true bright- 
ness distribution would show significant deviations, our results 
would be biased. Also the interpretation of the value found for 
the inclination (defined as the elongation of the Gaussians) then 
would be no longer straightforward. Unfortunately, such devia- 
tions can not be ruled out with the limited number of measure- 
ments at hand. However, Gaussian brightness distributions have 
proven to be a reasonable assumption for a wide variety of cen- 
trally heated objects. 
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Fig. 6. Histograms of the best fit parameters computed using 1000 Monte Carlo simulations for Haro 6-10 N (left) and Haro 6-10 S 
(right). 



5. Discussion 

In this section we compare our results with previous works and 
discuss which mechanism might have formed the binary system 
Haro 6-10. 



5.1. Large- and Small-Scale Structures in the Haro 6-10 
Environment 

The optical images of Haro 6-10 highlight some filamentary 
structures (as described in Section l4Tb which might have the fol- 
lowing explanations: 1) it is circumstellar material with an arc- 
like shape simply illuminated by the central star, or 2) it is the 
projection of the scattered light of a cavity in the circumbinary 
envelope, or 3) it represents the circumstellar material shocked 
by an outflow. 

The existence of a cavity in an en velope has been sugge sted 
around a similar system T Tau N, bv iMomose et al.l (1 19961) and 
IStapelfeldt et al.l (Il998ah . In this picture, a wide-angle wind cre- 
ates obUque shocks when interacting with the molecular en- 
vironment. On the other hand, the existence of an outflow is 
supported by the resolved extended emission of H2 at 6.5"in 
the south direction and i s consistent with the filament structure 
(iDoppmann et al.ll2008l) . iDevine et all (Il999l) already identified 
three outflow systems in their narrow-band optical images. In 
addition to the large-scale structure two separate outflows closer 
to Haro 6-10 were found. This suggests that both components 
are producin g outflows. Based on integral-fie ld spectroscopy and 
polarimetrv iMovsessian & MagakiarJ (Il999h suggest the north- 



ern source as the origin of a jet with a length of 5"superimposed 
on the bright reflection nebula. They propose that the disc around 
this component is seen edge-on, while the sout hern cornponen t 
is surrounded by a less inclined disc. Stoiimirovic et al.l (120071) 
report on a parsec-scale Herbig Haro flow with clumpy red- 
shifted gas pointing to the north-east. The blueshifted south- 
western part of the flow is contaminated by an unrelated fore- 
ground cloud. IBeck et a l. (2010) found a bulk of Bry emission 
from the two sources together with an extended, fainter emission 
from the southern component. The same extended emission was 
measured also in the [Fe 11] line and is thus clearly determining 
the direction of an outflow originating from this source. 

At 1.2 mm the emi ssion was again resolved with the PdBI 
(iGuilloteau et al.l [201 ih and the millimetre flux is almost the 
same for the two components. Under the assumption that the 
disc is optically thin at these wavelengths, the mass of the 
discs can be directly derived from the fluxes. The measured 
fluxes thus imply that the discs of Haro 6-10 N and S have al- 
mo st the same mass. The m easured large extinction associated 
by 'Stoiimirov ic et al.l (l2007h to a foreground cloud is also in 
agreement with the envelope we introduced to explain the mid- 
infrared absorption spectrum through both components of the 
binary system. Our findings on the circumstellar discs around 
the two sources of the binary system are consistent with the pic- 
ture drawn by the above studies. Also the misalignment of the 
discs and the suggested high inclination of the northern compo- 
nent is confirmed by our measurements. In particular, the almost 
face-on disc could explain why the emission of the spatially re- 
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solved Bry emission is observed blueshifted: the coiresponding 
redshifted emission might be covered by the face-on disc. Since 
the position angles of the discs can not be well constrained by 
our new data, we did our modelling (Section 14.3.2b with fixed 
position angles (T. Beck, private communication) as well as with 
the position angles as free parameters. Both approaches lead to 
the same results. 



5.2. Comparison to previous work 

The main results of our study are that Haro 6-10 is embedded in 
a common envelope (which causes the mid-infrared absorption 
feature in both components of the binary system), and that the 
two discs of the binary system are misaligned. 

Observational studies of the disc inclinations in young binary 
systems have been carried out in the last years, using different 
techniques: polarimetric and direct observations of protoplane- 
tary discs or indirect observations via the position of the jet from 
the protostellar object, which is expected to be perpendicular 
to the disc. Recently, also long-baseline infrared interferometry 
was used to measure the geometrical properties of circumstellar 
discs in binary systems. 

Polarimetric studi es on the d i sc ori e ntations ha v e bee n 
carried out by, e.g ., Monin et al.1 (Il998h : IWolfetal.1 (l200lh : 
iMonin et al.l (l2006l) . They all observed T Tauri sta rs and found 
that d iscs tend to be aligned in young binaries. Jensen et al. 
(l2004l) . in particular, found a tendency for binary system discs 
to be nearly (but not exactly) aligned with each other, and for 
those in triple and quadruple systems to be misaligned. 
There are also cases of misaligned discs found by obser- 
vations of misaligned jets from protostellar objects (e.g., 
Davis et al] 1994). infe rred jet precession (e.g., Eisloffel et al. 



1996 1: iDavis etai1ll997l). and d 
1998t IStapelfeldt et al."l998b') 



direct observations (e.g., iKoresko 



A recent interferometric study 
of T Tau showed that the discs in this triple system are mis- 
aligned ( Ratzka et alJ l2009): The circumstellar disc of the north- 
ern component is seen almost face-on, while the southern com- 
ponents are surrounded by higher inclined discs. Especially the 
disc around T Tau Sa, the most massive star in the system and 
a prototypical IRC, is seen almost edge-on. The formation pro- 
cess of T Tau thus has to have been highly dynamic. However, 
these cases of misaligned discs do not represent yet a statistical 
significant sample. 



lower However, disc-assisted capture is effici ent for dense envi- 
ronments and massive stars with large discs (iMoeckel & Ballvl 
2006, 2007b a). 

The most probable scenario is the formation via fragmen- 
tation of two different parts of the collapsing molecular cloud 
combined with other dynamical processes related to the cloud or 
the protostars. 

However, no final conclusion about the formation scenario 
can be drawn yet. More about the dynamics of the system has 
to be known. Nevertheless, our observations of Haro 6-10 show 
another system that has misaligned discs. Many of the binaries 
with infrared companions might harbour an edge-on and a face- 
on disc. 



6. Conclusions 

In this paper we presented a multi-wavelength study of the young 
binary system Haro 6-10 combining optical images obtained 
with HST/WFPC2 with the diffraction-limited images obtained 
with VLT/NACO in the near-infrared. In the mid-infrared we 
combine the spectroscopic and interferometric observations ob- 
tained with T1MM12 and VLTl/MlDl, respectively. This unique 
collection of data enables us to characterise the large-scale and 
small-scale structures of Haro 6-10 and to draw the following 
conclusions: 

1. Both components of the binary system Haro 6-10 are embed- 
ded in a common envelope. The envelope has a dust compo- 
sition similar to the interstellar medium. 

2. Each component of the system has a circumstellar disc-like 
structure typical of young stars. 

3. The discs are highly misaligned: the northern disc is seen al- 
most edge-on with respect to the line-of-sight and the south- 
ern component is surrounded by an almost face-on disc. 

4. The most probable scenario is the formation via fragmen- 
tation of two different parts of the collapsing molecular 
cloud combined with other dynamical processes related to 
the cloud or the protostars. 
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5.3. How Haro 6-10 formed? 

The misalignment of the discs in the Haro 6-10 binary system 
might suggest the formation mechanism of the system itself. 
Such a misalignment cannot be caused by an infall of external 
material with a different angular momentum onto the binary or- 
bit, because the a ccretion would have re-aligned the two discs 
(iBate et alJ l2000h . Only the infall on one of the components 
might lead to a tilt of that component, when the transported an- 
gular momentum is large enough. Also capture of the compan- 
ion is an unlikely explanation, because for solar-mass stars the 
capture rates are insufficient even when circumstellar discs in- 
crease the interaction cross section considerably (Helleri ll995h . 
For the central part of the Orion Nebula Cluste r (ONC ) with 
its stellar density of lO'* stars pc "* Bo ffin et al.l (Il998l) found 
a capture rate of only 0.1 Myr"', i.e. the probability to form a 
binary star within 1 Myr is only 10% . For environments like 
the Taurus-Auriga star-forming complex, with stellar densities of 
around 10 stars pc"^ the capture rate is two orders of magnitude 



Appendix A: Spectral Energy Distribution 

For completeness, we compile the spectral energy distributions 
(SEDs) of the two components of the binary system combining 
our new results in the optical, near- and mid-infrared, with the 
photometry from the literature. All the fluxes, together with the 
instruments, the beam sizes of the observations, and the refer- 
ences, are listed in Table I A. fl In Fig. lA.ll we show the SEDs of 
the individual components. 

In the wavelength range 12-1 000 jum the binary system is not 
resolved. But the level of the MIDI spectra reveal that most of 
the far-infrared flux originates in the northern component. 



Appendix B: Models 

In the paper the 'Two Black-Body Model' has been used to inter- 
pret the MIDI interferometric data. In the foUowing, we list the 
details of this model. 
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Fig. A.l. Spectral energy distribution (SED) of Haro 6-10. The empty diamonds points represent the photometric measurements of 
the unresolved binary. The fluxes of the single components are presented as red (Haro 6-10 N) and blue (Haro 6-10 S) filled dots. 
The continuous lines represent the MIDI spectra of both components. 

Table A.l. Spectral energy distribution with photometry from the literature (references in the last column) of Haro 6-10 N and Haro 
6-10 S. 



A 


Flux (N+S) 


Flux (S) 


Flux (N) 


Instrument 


Beam size 


Date 


Reference 


[jum] 


[miy] 


[mIy] 


[mJy] 




["] 






0.593 


0.438 ± 0.001 


0.437 ± 0.001 


0.001 ± 0.000 


HST 


0.5 


Nov. 1998 


this work 


0.792 


2.374 ± 0.004 


2.365 ± 0.004 


0.009 ± 0.000 


HST 


0.5 


Nov. 1998 


this work 


1.650 


39.2 + 3.4 


37.0 ± 3.4 


2.1 ±0.2 


NACO 


0.2 


Dec. 2003 


this work 


2.160 


166.3 + 12.6 


132.9 ± 12.2 


33.4 ±3.1 


NACO 


0.2 


Dec. 2003 


this work 


12 


16600 ± 2200 


- 


— 


IRAS 






Chandler etal. (1998) 


25 


37600 ± 3200 
58400 ± 9700 


- 


- 


IRAS 
IRAS 






Chandler etal.r 1998) 


60 


Chandler etal. (1998) 


100 


47000 ± 11000 


- 


- 


IRAS 






Chandler et al. (1998) 


350 


3420 ± 250 


- 


- 


JCMTAJKT14 


19 




Dent etal. (1998) 


443 


2370 ± 500 


- 


- 


JCMTAJKT14 


19 


Sep. 1991 


Chandler et al. (1998) 


790 


571 ±44 


- 


- 


JCMTAJKT14 


16 


Sep. 1991 


Chandler etal. (1998) 


1104 


179 + 17 


- 


- 


JCMTAJKT14 


19 


Sep. 1991 


Chandler etal. (1998) 


1260 


99 ±30 


- 


- 


JCMT/UKT14 


21 


Sep. 1991 


Chandler etal. (1998) 


1300 


90.5 ± 4.5 


46.7 ± 3.2 


43.8 ±3.1 


PdBI 


0.9x0.6 




Guilloteau et al. (2011) 


2700 


19.6 ± 1.0 


9.1 ±0.7 


10.5 ± 0.7 


PdBI 






Guilloteau et al. (2011) 


3400 


11.9± 1.7 


- 


- 


OVRO 


6 


1993/1994 


Hogerheiide et al. (1997) 


36000 


- 


1.10 


0.10 


VLA 


0.29x0.26 


Mai-ch 2002 


ReiDurth etal. (2004) 



The ID Gaussian intensity brightness distribution (I(a,/3), a 
andyS are the coordinates) is: 

1 /-4/n2p2\ 

I{a,/3) = -^==^ exp(— ^^J , 



yjnlAlnl® 



(B.l) 



where is the FWHM of the Gaussian, and p - -\Ja^ + jf- . 
The Fourier Transform of Eg. IB. H is: 

/ {Ti®rf\ 
^("'^)-xp(-i^), (B.2) 



where r - yiP^+v^ . 



The visibility of a Gaussian brightness distribution inclined 
by an angle / and rotated by an angle (which represents the 
position angle of the major axis) can be obtained by rotating the 
M v-coordinates (e.g., Berger &_Segransan 2007): 

(B.3) 
(B.4) 

(B.5) 



u — u ■ cos{ 



• + V ■ sind 



v' — —u- sin((p) + V ■ cos4 



\2 



■ cosiTp- . 

Here cos(iY represents a compression factor along the minor 
axis. 

The visibility of a 2D Gaussian brightness distribution is 
given by Eg. IB. 21 but with r defined according to Eg. 
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Fig. 7. MIDI spectrum, correlated fluxes and visibilities for Fig. 8. MIDI spectrum, correlated fluxes and visibilities for 
Haro 6-10 N. The red solid line represents the model best fit. Haro 6-10 S. The red solid line represents the model best fit. 



A circumstellar disc is approximated with two Gaussian 
components which emit as two black-bodies at different tem- 
peratures (Ti, T2)- Since the infrared spectrum is observed in ab- 
sorption, foreground dust characterised by its optical depth t is 
included. The black-body emitting at temperature T {B^{T)) and 



the corresponding total flux received (FJTJ) are: 



By(T) = ^expl 



(Jl:2_\ 

\k-T -\l 



Fa(T) 



n ■ q~ 

^2" 



Ba(T) , 



(B.6) 

(B.7) 
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where q is the radius of the emitting region (in At/) and d is the 
distance of the source (in pc). For each binary component, the 
mid-infrared spectrum is then computed as 



FUX)^\F,^i{T,)-e- 



+ FiJT2) ■ e 



■ cos(i) , (B.8) 



where F^^iiTi) and ^^,2(^2) are given by Eg. IB. 71 Since the disc 
is inclined the received to tal flux is multiplied by cos{i) (e.g., 
lBerper& Segransanll2007h . with / being the inclination of the 
disc. The function t{A) describes the extinction profile between 
8 and 13 |/m as measured in the direction of the galactic centre, 
and t(Ac) corresponds to the maximum of t{A). ti and T2 are 
the average extinctions computed in Section l4!2l for Haro 6-10 N 
and Haro 6-10 S, respectively. The visibility is computed as 



Vtot(/l)=/-Vgaussl+(l-/)-Vs 



^ ^auss2 



(B.9) 



where / - y^. ygaussi,2 are defined according to Equation ( IB.2b . 
with 012 - r\;ild, where ri_2 is the radius of the emitting region 
(mAU) and d is the distance of the source {mpc). Equation ( IB.9l l 
is valid only, if both F,i^{T\) and ^^,2(7^2) are axisymmetric to 
the same axis (one going through the star), which is assumed 
here. 

The correlated flux for each component of the binary system 
is: 

F,ou^ FUX)-VU^), (B.IO) 

where FiotiX) is given by Eq. IB. 71 and ytot(^) by Eg. IB. 91 
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